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|> ■ INTRODUCTION 
O 

o . 

(*C) . The normal state of the recently discovered high-T c superconductors [[!]], 

exhibits many exotic properties that are not yet fully understood. A picture 
\ that is gaining acceptance assumes that the normal state is not a normal 

Fermi liquid [2-5] but rather it shows unconventional properties somewhat 
similar to those of the Hubbard model in ID ||, namely: i) an unrenor- 
malizable Fermi-surface phase shift, which implies a vanishing wavefunction 
renormalization constant (Z), ii) a one-particle spectral density that shows a 
strong peak at the Fermi energy, hi) a Fermi surface that obeys Luttinger's 
^ ' theorem 0, and, iv) charge and spin separation. 

Only in very few cases has a quantitative check of these properties been 
^ undertaken [8-15]. In particular, Z has been calculated by exactly solving 

either the t-J |y| or the Hubbard Hamiltonian in 4 x 4 clusters. Charge 
and spin separation has also been recently investigated on finite clusters of 
the square lattice [[13|] . Unfortunately the behavior of these properties in the 
thermodynamic limit cannot be investigated at present by means of exact 
calculations. 



a 



In previous papers |l6l we have presented an Unrestricted Hartree-Fock 



(UHF) calculation of Z for the one-band Hubbard model in the square lat- 
tice. UHF allows to consider clusters large enough to investigate the scaling 
of Z with the system size. The results indicate that the system shows non- 
conventional behavior (Z=0) near half-filling. For a given value of U, Z 
vanishes at low dopings, and becomes finite as doping is increased. For very 
large U the value of Z remains finite for all dopings, excluding half-filling. In 
this work we discuss these results outlining the reasons why we think that 
UHF may provide trustable information on Z. We also show that the results 
are not changed if the wavefunction is written as a linear combination of 
UHF solutions each one centered on a cluster site (this procedure signifi- 
cantly improves the UHF results for the total energy). Finally we present a 
calculation of Z for the two-band Hubbard model; again, near half-filling, Z 
scales to zero with the system size. 

THE HUBBARD HAMILTONIAN AND THE UHF APPROXI- 
MATION 



The Hubbard Hamiltonian for one and two bands can be written in the 



general form Wl 



H = H + Hi, 



(1) 



h = J2 E i n ic - J2 kA c j<7, (2) 



<ij>cr 



H 1 = J2 Uin^-nu (3) 

i 

where the indexes i,j run over the atoms of the Cu02 planes in the two-band 
case, or over atoms at sites of the square lattice in the one-band case. The 
operator Cj CT destroys an electron with a z-component of the spin a at orbital 
i, and tij is the hopping matrix element between orbitals located at atoms i 
and j (the symbol < ij > denotes the sum over all pairs of nearest neighbors), 
Ei are the orbital energies and U{ is the intrasite Coulomb repulsion. 

The most general effective Hamiltonian, within the UHF approximation, 



can be written as |Tq, 19 , 



^ e// = E(44i)f < 7 > ^ <c ^>* 




Ui(< n iT >< nn > - < c^ Cii > ) (4) 



If spin flip terms characterized by a non-zero value of < cLc^ > are ig- 
nored, the standard Hartree-Fock approximation of the Hubbard interaction 
is recovered. In this work we have only considered solutions with a single 
non-zero component of the local magnetization for the following reasons: i) 
near half-filling the vortex solutions described in |18| , lie at much higher 
energies than the Ising solutions (magnetic polaron or extended), and, ii) at 
moderate dopings vortices become more competitive but only for an even 
number of particles. In general, solutions with transverse magnetization are 



THE WAVEFUNCTION RENORMALIZATION CONSTANT 



The wave-function renormalization constant is given by 

^ <M/(AQ|cA(iV+l)> 

< ^(N + l)\clc a \q(N + 1) > 

where a stands for the appropiate quantum numbers and the \l/'s are the 
ground state wave-functions. In calculating Z we proceed as follows. We 
start from the Slater determinant for N electrons and empty the one-electron 
state of highest energy. This is equivalent to identify the operator c£ in Eq. 
(2) with 



c 



1 -J2lUE n )cl, (6) 



where E n and a are the energy and spin of the empty UHF level of lowest 
energy, and 7k (T (-E'n) are coefficients obtained in the UHF calculation. For this 
choice of the operators in Eq. (5) the normalization constant in this equation 
is equal to one. Simultaneously, we introduce a small random distortion 
(maximum modulus of 10 -6 ) on the starting mean field, both on local charges 
and magnetizations. The addition of this distortion has two objectives: a) 
to break the degeneracy of the filled state of highest energy as occurs in 
cases such as half-filling, and, b) to draw the iteration process towards the 



state of lowest energy. Then we initiate the iteration process [18 , |19|1 until 
selfconsistency is achieved. 

The square root of the wave-function renormalization constant Z is given 
by the overlap between the resulting Slater determinant and the initial one, 
which can be calculated through the following expression, 



< tf(JV)|$(iV) >= 



(7) 



< ^jv|01 > ■ ■ • < 4l N \<pN > 

where *f>(N) and $(iV) are two Slater determinants for N particles, and ipi 
and (pi the corresponding monoelectronic wavefunctions. 

One-band Hamiltonian 

The one-band Hubbard Hamiltonian has been solved on finite clusters of 
size L x L (with L up to 16). For the parameters in Eqs. (1-3) we took, 
tij — 1, Ei — and Ui = U. The ground state for intermediate values of U 
is, in the case of a single hole, the magnetic polaron, whereas for many holes 
it corresponds to magnetic polarons disorderly distributed through the whole 



cluster. At U = oo the UHF ground state is ferromagnetic for all fillings [116 . 



Figure 1: UHF results for the wave-function renormalization constant Z 
vs the cluster size (L x L) for the one-band Hamiltonian and 6.25% doping. 



Our results for Z at half-filling clearly indicate that it slowly decreases 
with L for all U. The UHF results for the 4x4 cluster are in good agreement 
with the exact results reported in Refs. [ [T2|] . The scaling of Z to zero with 



the system size, can be understood as follows. The changes in the charge and 
spin configurations that result from emptying one of the degenerate levels at 
the top of the lower Hubbard band are non-uniformly localised along lines. 
Thus, as the one-hole selconsistent solutions are localized, it can be expected 
that the overlap between the two should decrease with the size of the cluster. 
It should be noted that as at half-filling the system is an antiferromagnetic 
insulator, strictly speaking it does not make sense to talk about its Fermi or 
non-Fermi like character. We note, however, that homogeneous HF solutions 
would have led to a finite Z even at half-filling. Thus, the present results 
clearly indicates that for a small, still finite, density of holes the system 
should show non-conventional behavior. 



Figure 2: Schematic phase diagram for the one-band Hubbard model ac- 
cording to the results for Z. The parameters are the fraction of holes referred 
to half-filling (x) and U. Z vanishes in a region close to half-filling. 

As regards finite doping, we have analyzed the case of 6.25% holes. The 
results of Figure 1 indicate that, for 6.25% of holes, Z scales to zero for U = 
20 whereas it remains finite for U = 10. This defines a transition from Fermi- 
to Luttinger-like behavior at an intermediate U. Thus, we may expect that 
at large enough U the constant Z will again be finite. At U — oo the UHF 
results show that Z remains finite at all dopings, except at half filling (this is 
a consequence of the ferromagnetic character of the UHF ground state away 
from half-filling). On the other hand, in the dilute limit the results of Ref. 
[ p0|| indicate that Z is finite for all U. 

The numerical data for Z, in the cases for which it decreases with L, can 



\[Z = aexp(-MnL) (8) 

This equation coincides with that suggested years ago by Anderson ^], PT[ |, 
and indicates that there exists an unrenormalizable phase shift at the Fermi 
energy. Actual results for b show appreciable errors due to the difficulty of 
the calculations. In all cases analysed here, b ~ 0.5. The resulting phase shift 
S, obtained by comparing Eq. (8) with that in refs. || |2lH , is 5 ~ or 



by introducing the Fermi wave number, 5 ~ 0.9kpa, in excellent agreement 
with the result of 0. 

The previous results suggest the phase diagram for the one-band Hubbard 
model depicted in Figure 2. The region where the system is expected to show 
non-conventional behavior is close to half filling and finite U. For a given 
finite U, the system will first (low doping) be a Luttinger liquid and become 
a Fermi liquid as doping is increased. 

Two-band Hamiltonian 

We have chosen the following set of parameters, E p — Ed = 4.0, Ud = 6.0, 
t p d = 1.0, and the remaining parameters equal to zero. This choice gives, 
for instance, a reasonable fit of the experimental data for the magnetization. 
Clusters containing up to 12 x 12 unit cells, have been considered. Both 
electron and hole dopings have been investigated. For a single hole the ground 



state is, in both cases, the magnetic polaron [19]. The results for Z at half- 
filling are shown in Figure 3. As in the case of the one-band Hamiltonian, Z 
decreases with L as in Eq. (8). The scaling of Z to zero is very similar for 
both types of doping. 



Figure 3: UHF results for the wave-function renormalization constant 
Z vs the cluster size (L x L) for the two-band Hamiltonian at half-filling. 
Results for both hole and electron doping are shown. 



BEYOND UHF 



The magnetic polarons, although show the appealing features of all mean 
field solutions Jl8|, [19|, have many drawbacks. The two most outstanding are 
the non-uniform distribution of the excess charge [13-15], in contrast with 
the exact solutions [5-9], and an energy, significantly higher than the exact 
one [I5|. In particular the former could be the origin of the behavior of Z 
discussed above. The usual way followed in quantum chemical calculations 
to go beyond HF is the so-called Configurations Interaction (CI). In the 
case of the magnetic polarons a particular CI seems rather clear, namely, the 
one that allows the localised polarons to move through the whole cluster. We 
have explored this CI for a single hole, by writting its wavefunction as a linear 
combination of Slater determinants for the single polarons (SP) localised on 
all cluster sites. The exact Hamiltonian has been solved in this basis set. We 
call this solution a multipolaron (MP). 

Thus, we write the wavefunction of a single hole (N s — 1 electrons) as 
a linear combination of single polaron Slater determinants centered on all 
cluster sites 



* MP (N a - 1) = E a i&( N s ~ !) ( 9 ) 

i 

where MP stands for multipolaron, and $ l (iV s — 1) is the wavefunction 
for the magnetic polaron centered on site i. The coefficients are obtained 
through diagonalization of the exact Hamiltonian. 

The ground state wavefunction for a single hole includes all the SP wave- 
functions with the same weight, and, therefore, gives an excess charge uni- 
formly distributed throughout the whole cluster, as in the exact case. Note, 
however, that the S z component of the single polaron (S z = 1/2) is preserved 
by the multipolaron approximation. The results for the energy necessary to 
create one hole (the difference between the energies of the system with one 
and zero holes) in the SP and MP approximations, are reported in Table I. 
As expected, the MP approximation lowers the energy of the single polaron. 
The most important changes are found for intermediate values of U. The 
correction rapidly decreases with U, and tends to zero for infinite U. On the 
other hand we note that the correction does only slightly depend on the size 
of the cluster. This can be easily understood by noting that the interaction 
between polarons rapidly decreases with the distance. 

We note that although the MP solution restores the translational symme- 
try of the Hubbard hamiltonian, features associated with the initial, localized 
solutions, remain; in particular the behavior of Z is not changed. This is so 
because, in our CI scheme, each of the solutions which enter in the final wave- 
function contains a "continuum" of levels which is substantially modified by 
the addition of an extra electron. The orthogonality catastrophe which is re- 
sponsible for the absence of coherent quasiparticles within the Hartree-Fock 
approximation is unchanged in the final wavefunction. 



CONCLUDING REMARKS 



Table 1: Energy (in units of t) required to create a hole in the single polaron 
(SP) and multipolaron (MP) approximations, for three cluster sizes (L) and 
several values of the intrasite Coulomb repulsion (U). Exact results are taken 
from Refs. [7,8]. The HF solution for the 4x4 cluster and U = 4 is not the 
single polaron (see text). 

LxL 





exact 




SP 






MP 




u 


4x4 


4x4 


8x8 


12 x 12 


4x4 


8x8 


12 x 12 


4 


-1.0434 


-0.7352 


-0.7317 


-0.7270 


-0.7942 


-0.9219 


-0.9366 


6 


-1.4148 


-0.8956 


-0.9031 


-0.9033 


-1.2143 


-1.2413 


-1.2419 


8 


-1.6784 


-1.0719 


-1.0793 


-1.0793 


-1.3731 


-1.3940 


-1.3941 


10 


-1.8640 


-1.2092 


-1.2190 


-1.2190 


-1.4662 


-1.4868 


-1.4868 


12 


-2.0015 


-1.3146 


-1.3267 


-1.3266 


-1.5326 


-1.5527 


-1.5527 


16 


-2.1954 


-1.4627 


-1.4767 


-1.4767 


-1.6249 


-1.6432 


-1.6432 


20 


-2.3281 


-1.5601 


-1.5741 


-1.5741 


-1.6871 


-1.7034 


-1.7034 



function renormalization constant for the one-band Hubbard Hamiltonian 
in the square lattice. Our results suggest that it may show non-conventional 
behavior in the thermodynamic limit, as discussed by previous authors [2-5]. 
The region of the parameter space (U and filling) where Z vanishes is nar- 
row and close to half-filling. The way in which Z scales with the size of the 
system makes this problem very difficult to investigate by means of exact 
calculations. 

It is worth mentioning that our route towards the results which show 
highly unconventional behavior is itself very conventional. We apply the most 
standard technique in the study of many-body systems: the Hartree-Fock ap- 
proximation. While the method, and the results, seem to be a contradiction 
of terms, we think that the opposite is true. Precisely because the method is 
so standard, it is biased in favor of normal (Fermi liquid) behavior. In fact, 
we recover such behavior in the low density limit, as we should. Also, the 
AF insulator at half filling is well described within our approximation. The 
main deviation of our results from conventional behavior is the proliferation 
of solutions which break translational symmetry, near half-filling. It is also 
in that region where we find Z=0. Restoring the translational symmetry by 
means of a CI scheme, while improving the results for the total energy, does 
not change the results for Z. 

Some results for the two-band Hamiltonian, indicate that this model may 
also show non-conventional behavior, no matter the type of doping. 
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